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ABSTRACT 


Stauffer (1973) developed a theoretical model for the prediction 
of air flow above ocean waves. In this study his results were extended 
with the addition of buoyance effects and two different lower boundary 
conditions. | 

Through numerical solutions it was possible to determine that: 
1) mean wind velocities at lower levels closely approximate the initial 
mean velocity profile for large wave numbers (k > 0.45m ~) or small 
wave amplitudes, 2) mean wind velocities at lower levels fluctuate more 
from the initial mean velocity profile as the wave number decreased or 
the wave amplitude increases, 3) solutions were very Sensitive to the 
level where the velocity gradient is computed in the lower boundary 
condition, 4) disturbance potential temperature at lower levels in- 
creases inversely with wave number, and 5) the presence of temperature 
fluctuations has an insignificant effect on the Reynolds stress or on 
the stream function. The numerical solutions show general agreement 
with the observational analyses of wave modified wind profiles by 
Davidson (1974). However, these comparisons are not concluSive because 


of the sensitivity in the lower boundary condition. 
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I. INTRODUCTION 


The forerunner of this study (Stauffer, 1973) used a 
theoretical numerical model to investigate wave-related 
momentum transfer. The model was used to obtain the Rey- 
nolds wave stress of the perturbation-velocity field in 
the atmospheric layer above ocean waves. The study was 
Similar to that of Yefimov (1970) and utilized an initial 
value numerical procedure similar to that of Newaee (Joe keys)) ae 
Equations were solved by utilizing central differencing in 
general but with the turbulence terms evaluated at the pre- 
vious time step. Solution separation was avoided by in- 
corporating the Matsuno (1966) finite difference scheme 
every 50 time steps. Experiments were run varying the in- 
dependent variables: wave number (k), dynamic velocity 
(U,), roughness parameter (20), and the coefficient of 
turbulent exchange (V). The measure of the error in the 
change of the stress in air (R) was computed at each time 
step and calculations continued until the stress reached 
meocoddyectace Condition defined by R < 1 x LO? Maximum 
Reynolds stress occurred at about one meter above the wave 
surface, and the largest momentum transfer was found to 
exist for the largest mean wind speeds, lowest wave numbers, 
and a turbulent coefficient of approximately 0.02m~/sec. 


The magnitude of the Reynolds stress increased as the 
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mean wind velocity approached that of the wave phase 
Speed in the lowest levels of the boundary layer. Thus 
when the wave velocity is close to the mean wind velocity 
in the lower levels, the disturbance energy is the 
greatest. 

The emphasis in this study is on the changes in the 
Mean wind which are caused by the’ wave stresses which 
were computed by Stauffer (1973). His basic model is 
modified to include buoyancy effects which arise during 
non-neutral conditions. Also the lower boundary condi- 
tion was changed to include the mean wind shear. Two 
approximations to this boundary condition were investi- 
gated. Experiments were carried out in which the follow- 
ing quantities were varied: wave number (k), grid size, 
convergence criteria, wave amplitude (a), and lower 
boundary condition. It was found that the mean wind de- 
parted from the logarithmic profile in the lower region 
when the wave number became small and when the wave ampli- 
tude was large. Some of these variations were consistent 


with the observations of Davidson (1974). 


iy 





Il. DEVELOPMENT OF THE MODEL 


Development of the model is Similar to that of 
Stauffer (1973). The Boussinesgq equations for two dimen- 


Sional motion take the following form: 














Eee WV CNN Soom ee gt Fe (2.1) 

se + Ve VO = - Z_ (u'O') - o (w'@") and (2.2) 

oo + eM = Oa, (2.3) 
where 

st o_ (-u'?) +r 2 (-u'w') , 

Fi — (=uaw) + s- (-w'7) ’ 


Here T(p /p)~ is the total potential temperature, ox is 
the mean potential temperature for the whole region, R 

is the universal gas constant, oes 1s the specific heat 

at constant pressure, x is the horizontal axis, and z is 
the vertical axis. Equation (2.1) is the equation of 
motion; equation (2.2) is the first law of thermodynamics, 


and equation (2.3) is the continuity equation. The 
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quantities u', 


‘ 
Woes 


and 


6! 


are turbulent fluctuations. A 


straight line above a variable signifies the stationary 


part (time averaged) of each variable, 


and a wavy line 


represents the non-stationary part (space averaged). The 
velocity (W) and the potential temperature (9) are 
broken up as follows: 
W = (U(x,t) + u(x,z,t)) |i + wh , ana (onan 
Or = fomz,t) + 8 (x,z,t) . (225) 


Following Yefimov 


mated as follows: 





-u'w' = T= K 
“wO" = Ky 55 
=e Gg? = Ses 


Here K and Kp 


for the mean fields, 


(1970), 


the 


the wave-caused fluctuations. 


The following vorticity equation is 


| 5 e Vx with equation (2.1): 


oF o€ 

ve + U ae + w 
2 

02 





3 


turbulent fluxes are approxi- 


and V and V 


(2.6) 


(2.7) 


(2.8) 


(229) 


are the coefficients of turbulent exchange 


are the coefficients of 


6 


formed by taking 





where 


_ du OW 
gegt- 2. | (25109 


in this derivation, equations y(2.3)), (2.4), (2.5), (2.6), 
and (2.7) have been used. When equations (2.4), (2.5), 
(2.8), and (2.9) are introduced into equation (2.2) it 


becomes: 








30 36 a0 30 90 9 30 06, 38 
Bim 6 Coe ee fae = Fe nas! ' 6 im © Wap) 
ox dz 
(23.22) 


These derivations have used the following assumptions: 


0 (2523) 


Take the x-average of equations (2.10) and (2.12): 


Zz 
Bre a aU 
z 
90 90 _ 2 30 
“ox * “az ~ 32 ‘Xe a2) ° ahaa 


Subtract Kzema),) and (2.15) from (2.10) and (2.12), respec= 
tively, and neglecting products of the fluctuation quanti- 


ties gives: 














2 2 
3 3 

d& a& dU re) du aw du od W g 3008 
eta + W =— _(vi=x— + =) I] -vli + —~—]J- = =~, 
dot ox \2 5 2 dz ) 54232 542 00 ax 

(2266) 
96 36 90 _ 0°68 aa 
7t ste Sone ote ae ver 5 + ; 5) e (2.197) 
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Since the motion is non-divergent (equation (2.3)), a 


stream function can be defined such that: 


oe ye Ob 


a o5 a as dx e (2.18) 


When these relations are introduced into (2.16) and (2.17), 


they become: 


2 2 3 
aay, 2 Ly BY 

















a7 w ow o7u 
+ ]+< 
Je An” ae Reeae \x> Ox ee 
2 2 2 4 4 : 
‘) dp dW dW dp g 00 
Brea, LV emmy) by ee ey ee (2.19) 
x 
Oz OZ ox Ox OZ ox O 
99 96 , dw 20 a29 . 926 
Bt wok  sscnae Oe oe zi Beuiey 


In this study the motions are forced by a surface wave of 
wave number (k) and phase speed (c). Therefore, the stream 
function and the potential temperature departure are 


written as: 


|} 
ll 


A(z,t) cos [k(x-ct)] + B(z,t) sin{k(x-ct)] ,(2.21) 


@D 
lI 


D(z,t) cos [k(x-ct)] + E(z,t) sin{k(x-ct)] .(2.22) 


Substitute relations (2.21) and (2.22) into (2.19) and 
separate the cosine and sine terms. Equating the coeffi- 


cients of the cosine terms gives: 

















2 Z 2 
ee a= kc) (Bk — ek 
at 2 2 
OZ OZ OZ 
2 2 2 
Moe a ype veo [eee + AK?) - oo (2.23) 
O. 5 2 0 
oz OZ OZ fe) 
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The coefficients of the sine terms give: 

















2 2 2 
328 - x7B) = -k(u-c) (ak?- 23)- ax 29 
OZ oz OZ 
2 2 2 
+ 2 v2 ® + Bx’) ES Pgs 
2 Q 
dz OZ oz re 


If the same procedure is carried out with equation 


the cosine coefficients give: 


dD _ 30 Pod 
ne -k(U-c)E + 37 «KB = Vp [Dk ea 


The sine terms give: 


~ 


0E _ a0 2 
ae k(c-U)D + a> kA - vo lEk = 





Substitute (2.21) into (2.14) yields: 


[A 








ie ee Ae dz 


Oz 


Nin 


(2.24) 


C220 )0, 


(2525) 


(2226) 


(2227) 


Equation (2.15) will not be used in this study because it 


turns out that the temperature effects are very small. 


~ 


Therefore 0 is equal to its initial value in this study. 


The equations (2.23) through (2.27) form a complete set 


in the variables A, B, D, E, and U. 
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gE BOUNDARY CONDITIONS AND INITIAL CONDITIONS 


Yefimov (1970) defined the boundary conditions for the 


lower boundary of the turbulent layer as: 


U(N,t) + (xo “akc cos [k(x-ct)] , (3-1) 


akc sin [k(x-ct)] , (3-2) 


w (x,n,t) 
where the height of the wave surface is given by: 


Nn = a cos [{k(x-ct)] . “(3-23) 


The deep water wave phase speed is: c = Re ygel » (3.4) 


If (3.1) and (3.2) are in terms of the stream function 


(2.18), they can be written as: 


) 
Uti, t) + ot (x,nN,t) = ake cos [k(x-ct)] , (3.25) 
ow ; 
ies (x,N,t) = ake sin [k(x-ct)] . (3.6) 
The stream function and its derivatives at Zz = NQ can be 


approximated by the values at z= 0, but U(n,t) cannot be 
approximated by U(0,t) since U has a very rapid change 
near the surface. Therefore U(n,t) is expressed with a 


Taylor series expansion as: 


_ aU 
U(n,t) = U(O,t) + a Z (33:7) 


Substitute (2.21) into (3.5) and use (3.7) and (3.3) which 


gives: 
oh WO, coe SenUee ee ~) ; (3.8) 
SB (0,t) anes (3.9) 


Ay 





The condition U(0,t) ~ O has been used. When (2.21) is 


Substituted into (3.6), the following are obtained: 


A(0O0,t) = ac , (32060 ) 


i 
© 
9 


B(0,t) (sei) 


The following boundary condition is used for the potential 


temperature fluctuation: 
6(x,0,t) = 0 (Se ) 


This condition isS not as accurate as the velocity conditions 
above, but the results shown later indicate that the tem- 
perature effects will be small in any case. Substituting 


We22) anto (3.12) yields: 
D(O,t) = E(0,t) = 0 (3473) 
In order to close the problem, the following conditions are 


imposed at the top of the region: 


Gat) = Su, t)= B(H,t)= = (H,t)= MCE ee (Hye 0 95514) 


where z = H is the upper boundary. 


The initial conditions are given by: 


A(z,0) = ac sinh[k (H,-z)]/sinh KH, (SSS) 
Boz, 0) c= 0 , (3.16) 
Dtz,0)) = 10 ,; (3.17) 
E(z,0) = 0 , (34 Lo) 
Uz . 
and U(z,0) = amma a) oman 1 ee (3.19) 
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where: Hy = H-Az, U, = dynamic velocity, Zee roughness 
parameter, K = 0.35 = von-Karmen's constant, and a = wave 
amplitude. The equations give an ferotati onal initial 
State with no temperature fluctuations. The initial mean 
wind is logarithmic. 


The following fields which are independent of time 


take the following form: 


K = K U,z , (35200) 
~ Tis 

vA 
QO = zm Uin Pe 1) , . | aeeeee) 


where T, is the temperature scale. 
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IV. NUMERICAL PROCEDURES 


Centered finite differences are introduced into equa- 
tions (2.23) through (2.26) with the exception that the 
turbulence terms are evaluated at the previous time step. 
The left hand side of the equations are solved for the 
time tendencies with the Guass elimination technique which 
is described in Richtmyer (1957). Centered time differenc- 
ing is used except that the solution is restarted epee 50 
time steps with the finite difference scheme developed by 
Matsuno (1966). This procedure eliminates solution 
separation. 

Note that the value of U given in equation (3.19) with 
the coefficient given by (3.20) leads to a zero value for 
the right hand side of (2.27). Therefore U can be expressed 
as follows: 

U 


: = u 4.1 
— 7 aed 1) + Ww(z, t) e (4.1) 


it 


U(z,t) 


Substitute (4.1) into (2.27) and integrate three times 


with respect to z which yields: 


een g fra, (4.2) 


where T = -uw = - =k (——A - ~B) (4.3) 
OZ OZ 


is the Reynolds wave stress. In this derivation the right 
hand side of (2.27) was rewritten and the coefficient K was 


replaced by the constant Vv. 
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The boundary conditions (3.8), (3.9), (3.13), and 
(3.14) are applied ina straightforward manner with a 
one-sided derivative when required. 

The time integration of the system of equations is 
accomplished by first solving the finite difference form 
of (2.23) through (2.26), and then computing the new U 
from aaaecon (4.1). Integration:‘is continued until a 
Steady state is reached. The following quantity is used 
as a measure of the steadiness of the solution: 


= 
2 


veal 
7 Az 
ae 


i 
The integration is continued until Ra > R, where Ri will 


take on three different values. 


Pad 





Ve. NUMERICAL SOLUTIONS 


The following constants were used in all numerical in- 


tegrations: 


Cc 
Il 


0.167 m sec” + ’ 


Zo = 7.9 x ae Mass C5 sls) 


Ow 24 ie Sica e 


Hl 


Three forms of the boundary condition (3.8) were used 
in these numerical integrations. Stauffer (1973) used the 


following condition which neglects the ee term: 


Oz 
S30, t) = -ake , (5.2) 


The direct use of (3.19) in equation (3.8) leads toa 
Singularity when applied at z = 0. Ina private communica- 
tion, Professor R. Davis has suggested that this problem 
can be alleviated by computing the derivative at a false 


level z = z In this case the boundary condition becomes: 


ln 


U, 


oA =o ee 
a, (Ort) = -al[kec + ce a Zz) ° C33) 


This formula reduces to (5.2) when z, is very large. This 


i 
procedure does not take into account changes in the mean 
wind near the surface through u. 

A third boundary condition was used which changes 


. 


as the mean wind changes. This takes the form: 


Iso 





OA U(Az,t) - U(0,t) 
yp (Ort) afkc + ie | (524) 


Numerical solutions were obtained varying six indepen- 
dent parameters as shown in Table I. These were wave 
number (k), wave amplitude (a), measure of the error in 
the change of the stress in air (R_), false level (z,), tem- 
perature scale (T,), and the boundary condition PETG, U)- 
Experiments were carried out as exhibited in Table II for 
the different values shown in Table I. 

Figures 1 and 2 show the distribution of mean wind 
velocity (U) with height (z) for series 1 and 2, respective- 
ly. These experiments use the lower boundary condition 
used by Stauffer (1973). Both figures show that the mean 
wind is close to the logarithmic profile for the largest 
wave numbers. They also show large variations as the wave 
number decreases. These departures are principally in the 
lower layers below approximately six meters. The smallest 
wave number (k = 0.050 m7 +) shows a much larger velocity 
than the logarithmic profile with a first maximum at approxi- 
mately 1.5 meters. The character of the solutions for the 
two values of Ra are similar, but for particular wave 
numbers there are considerable differences. 

Figures 3 and 4 represent series 3 for wave numbers 
k =,0.35 a and k = 0.20 m+, respectively. These actually 
represent the Reynolds wave stress (T) as a function of 
height. Curve six in each figure is for the case Zz, = 


1 


infinity which corresponds to the boundary condition used 
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Figure l. 
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Distribution of the mean wind velocity (U in 


with height 


DO fase. 2) 


(z in meters) 


Oped Seine) 


utilizing series 1] for 


respectively to: 1) initial mean 
0.375 m-l; 4) 0.275 m1 


5) 0.175 m~~; 6) 0.100 m=!; 7) 0.050 m1. 
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Figure 2. Distribution of mean wind velocity (U inm see os 
with height (z in meters) utilizing series 2 for wave number 
Gy) equaw™ respectively, to: l)inittalemean veloecuty pro- 
file; 2) 0.475 m71; 3) 0.375 m71; 4) 0.275 m+; 5) 0.175m °; 
6) 0.100 mt; 7) 0.050 ml, 
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Figure 3. Distribution of the wave-caused stress (T in 
0...01 m2sec 2) with height for wave number (k) equal to 
0.350 m7t, and utilizing series 3. Variation in Z4 is 
Sway respectively, to: 1) @27025 m: 2) 0.675 m;: 

Spee ee 25 me 4) 0.075 m; 5) OF025 m: 6) infinity. 
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Figure 4. Distribution of the wave-caused stress (T in 
0.09 m2sec 72) with height for wave number (k) equal to 
O22010 m1, and utilizing series 3. Variation in 24 is 
equidieerrescrPeetivery=ro;: Py 27025 mm; 2) )0.675 me; 3) O.225m; ° 
‘eer oun eo O.O25 ms 6) infinity. 
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by Stauffer (1973). These two curves are very similar to 
the curves for zy = 2.025 m. As Z4 decreases, the curve 
changes rapidly in such a way that the maximum stress for 
Zz, = 0.025 m is at least one order of magnitude greater 


than the case 2, = infinity. These experiments show the 
sensitivity of the solutions to this lower boundary condi- 
ELOn . 

Series 4, 5, and 6, ween ise the finite difference 
estimate of the wind shear in the boundary condition, were 
carried out with three values of Rae The solutions for 
Ra = Le? and Ra = ome showed some significant differ- 
ences while the solutions for Ro = ome were nearly identi- 
cal roe those of Ra = Om Figures 5 and 6 show the mean 
wind profile for series 5 which employs Ro = oa These 
figures show that for k >0.450 a the wind profile is 
very close to the initial mean velocity profile. Also 
note the large departures from the logarithmic profile at 
low levels for k = 0.1m : or less. A disturbing feature 
of the solutions is the oscillation of the mean wind be- 
tween adjacent values of the wave numbers. 

Figure 7 shows the solutions for the mean wind for 
series seven. These solutions are the same as those of 
series 4 except that the grid size has been cut in half 
and the time step has been reduced to one-fourth. These 
solutions are generally similar to those in figures 5 and 
6, except that in the lower levels the wind is larger than 


the logarithmic value for k <0.42535 m These results 
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Figure 5. Distribution of the mean wind velocity (U in 
m sec7~l) with herghe (Z in meters) utilizing series 5 
fer wave number (k) equal, respectively, to: 1) initial 
mean velocity profile; 2) 1.000 m~1; 3) 0.500 m~1, 

4) 0.450 m7; 5) 0.400 m1; 6) 0.350 m1, 7) 0.300 m7}, 
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Figure 6, 
m sec7l) 

wave number 
velocity profile; 
0.100 mt; 


3) 


with height 





Distribution of the mean wind velocity 
utilizing series 5 for 
initial mean 
Oo, 150 m~1; 


(z in meters) 
respectively, 
0.250 m7l; 
0.050 ml. 
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Figure 7. Distribution of the mean wind velocity (U in 
m sec7~l) with height (z in meters) utilizing series 7 
for wave number (k) equal, respectively, to: 1) initial 
mean velocity profile; 2) 1.000 m1; 3) 0.42535 m-1; 

4) 0.400 m71; 5) 0.200 m+; 6) 0.175 m-1; 7) 0.150 m-1, 
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are qualitatively consistent with the empirical relation 
developed by Davidson (1974). His data show departures 
from the logarithmic wind profile which are proportional 


to 0.16 (= - 26.3). For this condition in these experi- 


* 
ments, this quantity is positive for all wave numbers 
shown except for k = 1.0m : for which it is slightly nega- 
Hive. This behavior is consistent with the curves given 
in figure 7. One further test was made in this series 
for k = 0.200 a but with Ro = 10m. In that experiment 
the maximum wind velocity in the lower levels was very much 
larger than any shown in figure 7. This further emphasizes 
the sensitivity of the solutions to the lower boundary 
condition. 

Figures 8 through 1l are comparisons of data already 
presented for certain wave numbers. Figure 8 shows identi- 
cal solutions for the reduced grid size for k = 1.0 =a 
These solutions have a slightly See velocity in the 
lower layers than for the logarithmic profile. Figures 9, 
21), and 11 which are for wave numbers k = 0.40 ae 
k = 0.20 ae and k = 0.15 a respectively, show ir- 
regular variations about the logarithmic profile for dif- 
ferent values of Rv but systematically higher values of 
wind speed for Az = 0.125 m. 

_ Figures 12 and 13 contain the solutions from series 8 
for wave numbers k = 0.40 m and k = 0.15 m7, respective- 


ly. In all other experiments the wave amplitude is given 


by: 
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Figure 8. Distribution of the mean wind Wekoecuty (Y in 
m sec”) with height (z in meters) for wave number (k) 
equal to 1.000 m7+1t, Comparison is made where: 1) initial 


mean velocity profile; 2) series 4; 3) series 5 and 6; 
4) series 7. 
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Figure 9. Distribution of the mean wind velocity (U in 

m sec. ~) with height (z in meters) for wave number (k) 
equal to 0.0400 m~1. Comparison is made where: 1) initial 
meaievelocity profile; 2) series 4; 3) series 5 and 6; 

4) series 7. 


36 





024 


O21 


018 


os 


003 





Figure 10. Distribution of the mean wind velocity (U in 
m sec7l) with height (z in meters) for wave number (k) 
equal to 0.200 ml, Comparison 1S made where: 1) initial 
Mean velocity profile; 2) series 4; 3) series 5 and 6; 

4) series 7. 
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ees tae Distribution of mean wind Veloecrty (U in 
m sec71) with hei ht (Zz in meters) for wave number (k) 
equal to 0.150 m ~~. Comparison is made where: ye tne al 


mean velocity profile; 2) series 4; 3) series 5 and 6; 
4) series 7. 
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Figure 12. Distribution of the mean wind velocity (U in 
m sec +) with height (z in meters) utilizing series 8 for 
wave number (k) equal to 0.400 m7. variation in wave 
amplitude (a in meters) is equal, respectively,to: 

1) initial mean velocity profile; 2) 2.00 m; 3) 1.00 n; 
aymeOrm>O8m; 5) 0. 25%m=!i .1/kK; 6) 0.10m. 
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Figure 13. Distribution of mean wind velocity (U in 

m sec~l) with Mereie «2 inemeters) "WwEldldzing sierdies 8 for 
wave number (k) equal to 0.150 m~l, variation in wave 
amplitude (a in meters) is equal, respectively, to: 
Meter al Mean velocity profile; 2) 2.00 m; 3) 1.00 m; 
memo mM = 2.t/ks 5) O.50 m; 6) 0.10 m. 
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a = 0.1 k (2.0) 


In these figures, the amplitude ranges from 0.10 m to 
2.00 m. The smaller values of amplitude gives wind pro- 
files which are very close to the logarithmic profile. 
The larger wave amplitudes lead to increasingly larger 
deeartures from the logarithmic profile. This is to be 
expected since u is proportional to the square of the 
disturbance amplitude (equation 4.2). 

Figure 14 contains the disturbance potential tempera- 
ture as a function of z, from series 9 with three values 
of k. These curves show that the largest temperature 
amplitudes occur at about one meter in height; above these 
Maxima the temperature decreases rapidly. The largest 
amplitude occurs for the smallest wave number. The 
numerical solutions showed that the stream function ampli- 
tudes were not significantly affected by the presence of 
temperature fluctuations. AE jiibisy ige tee be ilove sehe pu UN 
fOelomand To. — —-O.15. This behavior may be interpreted 
as follows. If it is assumed that the horizontal advec- 
tion is balanced by the vertical advection, then equation 


(2.17) becomes: 


eee 
U-c Oz wW ® CSG.) 





This expression is used to estimate the last term in (2.19). 
It can be seen that the buoyance term will be much smaller 


than the vorticity advection term with the possible 
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Figure 14. Distribution of disturbance potential tempera~ 
ture (0 in °K) with height (z in meters) utilizing series 
9 for wave number (k) equal, respectively, to: 1) 0.275 m7; 
2) 0.125 m-1; 3) 0.050 m1, 
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exception of the point where U = c. The numerical solu- 
tions show that the presence of such points does not change 


the conclusion that the buoyancy term is that unimportant. 
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VI. CONCLUSIONS 


In this study the numerical results of Stauffer (1973) 
were extended with the addition of buoyancy effects and 
two different lower boundary conditions. The buoyancy 
effects have no appreciable affect on the stream function 
and wind stress fields; however, the presence of buoyancy 
effects could affect the Reynolds wave stress through modi- 
fications in the mean wind profile. These effects were 
not considered in this study. 

In this thesis, emphasis was placed on modifications 
in the mean profile which arose from the Reynolds wave 
stress. The two new boundary conditions which were inves- 
tigated involved two methods of estimating a in one of 
the lower boundary conditions. This term was entirely 
neglected in the study by Stauffer (1973). In the first 
method the derivative was evaluated from the logarithmic 
wind profile evaluated at a false level (21). This pro- 
cedure leads to very large wave stresses near the surface 


when a became small. 





ou was evaluated by finite dif- 


ferences. Many solutions were performed with this 


In the second method, 


boundary condition involving different values of the wave 
number, convergence criteria, and the grid size. Those 
solutions showed that the wind profile was near the loga- 


rithmic profile for the larger wave numbers of order one. 
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Smaller wave numbers showed considerable fluctuation 
about the logarithmic profile in the lower layers with 

the tendency toward wind speeds larger than those given 
by the logarithmic profile. Figure 7, which employs the 
smallest grid size, shows this tendency clearly. However, 
further investigations with this profile suggest larger 
fluctuations if a finer convergence criteria is used, 

In general the computed wind profiles are consistent with 
observational data analyzed by Davidson (1974). 

Further investigations of this type must carefully 
consider the lower boundary condition which was shown to 
be very sensitive in this study. Perhaps a full two-dimen- 
sional treatment of the air flow over the wave surface is 
required. Also it is very important to obtain accurate 
expressions for the turbulent exchange coefficients at 
different. locations near the wave surfaces. These studies 
should be coordinated with careful analysis of observa- 


tional data. 
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